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Abstract For regulated exocytosis, we model the morphology and dynamics of the 
making of the fusion pore or porosome as a cup-shaped lipoprotein structure (a dim- 
ple or pit) on the cytosol side of the plasma membrane. We describe the forming of 
the dimple by a free boundary problem. We discuss the various forces acting and an- 
alyze the magnetic character of the wandering electromagnetic field wave produced 
by intracellular spatially distributed pulsating (and well observed) release and bind- 
ing of Ca 2+ ions anteceding the bilayer membrane vesicle fusion of exocytosis. 
Our approach explains the energy efficiency of the dimple forming prior to hemifu- 
sion and fusion pore and the observed flickering in secretion. It provides a frame to 
relate characteristic time length of exocytosis to the frequency, amplitude and direc- 
tion of propagation of the underlying electromagnetic field wave. 
We sketch a comprehensive experimental program to verify - or falsify - our math- 
ematical and physical assumptions and conclusions where conclusive evidence still 
is missing for pancreatic j3 -cells. 

Key words: Calcium oscillations; dimple formation; free boundary problems; fu- 
sion pore; Lorentz force; Maxwell equations; pancreatic beta-cell; plasma mem- 
brane; regulated exocytosis. 



1 Introduction 

1.1 On Our Heuristic (Suggestive) Use of Mathematical Modeling 

This chapter adds a few electromagnetic facts and mathematical theorems to the 
toolbox approaching the process of bilayer membrane vesicle fusion. We address the 
related geometric and dynamic aspects of the endocytotic-exocytotic cycle which is 
at the core of various discharge (e.g., secretion) and ingestion (e.g., drug intake) 
processes in animal cells. 

We begin with a caveat. From mathematical physics, quantum chemistry, and 
various fields of engineering design we are accustomed to perfect reliability of the- 
oretical calculations due to full understanding of the governing laws and full prac- 
tical control of the calculated processes. We have learnt, sometimes the hard way, 
from physics history that, when in doubt, we had better trust the theory and carefully 
designed elaborate experiments than first views and ad-hoc explanations. 

Clearly, the situation is different in theoretical biology and medicine. There, it 
seems to us, the main use of mathematical modeling is either falsification or ex- 
trapolation. By falsification, we mean the use of simple arithmetic or other more 
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advanced mathematical means to check and falsify common belief (like Harvey's 
mathematical microscope, see Ottesen [1|, or the harmonic analysis of Ca 2+ os- 
cillations in Fridlyand and Philipson 1201 ). By extrapolation we mean the packing 
of established phenomenology into a precise, intentionally simplified mathemati- 
cal frame work admitting series of computer simulations or analytical estimates to 
investigate the role of selected parameters (like the Silicon Cell, see Westerhoff et 
al. [2 1, the mesoscopic simulation of membrane-associated processes, see Shillcock 
O), or the compartment models for different pools of insulin granules in exocytosis 
preparations, see Toffolo, Pedersen, and Cobelli (4)). 

In this chapter, our approach is different. We simply ask: (1) Could it be that a 
highly localizable phenomenon like the lipid bilayer fusion of regulated exocyto- 
sis on a characteristic length scale of tens of nanometers has essential cell-global 
aspects on a characteristic length scale of hundreds and thousands of nanometers? 
(2) Could it be that the observed changes of the electrostatic plasma membrane po- 
tential accompanying regulated exocytosis and the corresponding Ca 2+ oscillations 
have an electro-magnetic character which requires a field-theoretic (Maxwell) ap- 
proach to the secretion process? 

We have good reason for our, at present still speculative but hopefully suggestive 
approach, both in the re-interpretation of more or less well-observed phenomena 
and in focusing on aspects which seem to us not sufficiently supported by common 
explanations. This will be explained below. 

We shall emphasize that the correctness of the electro-dynamical and mathe- 
matical modeling parts of the findings of this chapter depend on future experimental 
testing and biological validation. At the end of this report, the reader can find a com- 
prehensive list of experiments that really need to be done to confirm the relevance 
of our physical equations and mathematical modeling. Therefore, our models do not 
aim for instant clarification but rather set a scene for alternative considerations and 
future observations. That's what we understand by the heuristic (suggestive) use of 
mathematics. 



1.2 Electromagnetic Free Boundary Route to Fusion Pore Making 

For regulated exocytosis, we model the morphology and dynamics of the making 
of the fusion pore or porosome as a cup-shaped lipoprotein structure (a dimple or 
pit) on the cytosol side of the plasma membrane. One ingredient to our model is a 
free boundary problem for the dimple under the action of electromagnetic forces, in 
particular the Lorentz force acting on charged molecules of the cell's plasma mem- 
brane with decreasing capacitive reactances while forming the dimple. The force 
comes from a wandering electromagnetic field produced by intracellular spatially 
distributed pulsating (and well observed) release and binding of Ca 2+ ions. 

Our approach is based on variational principles and emphasizes regularity and 
singularity under the deformation process of the membranes. It explains the energy 
efficiency of the observed dimple forming prior to hemifusion and fusion pore and 
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the observed flickering in secretion. It provides a frame to relate characteristic time 
length of exocytosis (ranging between milliseconds in nerve cells and seconds in 
j3 -cells) to the frequency, amplitude and direction of propagation of the underlying 
electromagnetic field wave. 

We shall not address all the machines (both protein machines and lipid assem- 
blies) working together in making the structure and the composition of the fusion 
pore. Admittedly, conclusive evidence is still lacking of the critical character of the 
here described electromagnetic field wave for the well-functioning of the regulated 
exocytosis in healthy cells and the lack of secretion robustness in stressed cells. In 
germ, however, the present electromagnetic free boundary model gives various hints 
to future calculations, estimates, and in vivo, in vitro and in silico (i.e., numerical 
simulation) experiments. 



1.3 Plan of the Chapter 

In Sect. [2] we summarize several mathematical, electrodynamical and cell physio- 
logical facts which seemingly have been overlooked or discarded in the literature, 
but may in our perception add essential ingredients to a comprehensive understand- 
ing of the short very first phase of regulated exocytosis. In Sect. [3] we describe our 
model and the corresponding differential equations, force balances and cost func- 
tional In Sect.[4]we discuss regularity and singularity results. In Sect.[5]we present 
our preliminary conclusion; some hints regarding the question of what controls the 
speed of the process; and a review of experimental tasks and capabilities to test our 
hypotheses. 



2 Synopsis of Established Facts 

We describe the challenge of making the fusion pore; distinguish different math- 
ematical modeling approaches; and elaborate electromagnetic and geometric phe- 
nomena of the very first phase of regulated exocytosis, namely Ca 2+ oscillations, the 
corresponding slow and low frequent electromagnetic field wave, and the forming 
of a dimple in the plasma membrane prior to the bilayer membrane vesicle fusion. 



2.1 Membrane Fusion and the Fusion Pore Challenge 

In animal cells, membrane fusion between the plasma membrane and transport vesi- 
cles is fundamental for the secretion of macromolecules. In contrast, the opposite 
event, i.e. the forming of vesicles or endosomes from the plasma membrane is nec- 
essary for the uptake of macromolecules and nutrients. The latter process is known 
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as endocytosis, in which, in the suggestive words of a renown textbook "localized 
regions of the plasma membrane invaginate and pinch off to form endocytotic vesi- 
cles" (Alberts et al., [5]). The process of discharge of material, after collecting it in 
transport vesicles, is called exocytosis and is our subject. 

It happens by docking of the vesicle to the cell membrane through activity of 
several membrane-associated proteins, followed by vesicle membrane hemi-fusion 
and the making of a fusion pore in the membrane through which the material can 
be released to the exterior, see Fig.[T] As a rule, the two processes, i.e. endocytosis 
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Fig. 1 Schematic view of the bilayer membrane fusion event, from Koch et al. 1 6 1 (after Lentz 1 7 1), 
reproduced with permission of World Scientific Publishers and the authors. Shown are the basic 
molecular structures involved in the process of fusion pore making and membrane fusion. Note 
the lipid bilayers of the plasma membrane and the vesicle membrane, with the negative charges 
representing the charged heads of phosphatidyl serine. For clarity, only the transmembrane proteins 
syntaxin and synaptobrevin with the three aromatic amino acids on the luminal side of the vesicle 
membrane are depicted. Also included is synaptogamin, which binds Ca 2+ and is associated to 
both syntaxin and synaptobrevin. 



and exocytosis, seem to be balanced, the one cutting membrane pieces out of the 



6 



D. Apushkinskaya et al. 



cell membrane, the other inserting pieces. In both processes, the remarkable is the 
opening of the cell without pinching a hole. 

Typical examples of endocytosis are the intake of nutrition, signal molecules, 
viruses or drugs; typical examples of exocytosis are secretion of macromolecules 
such as hormones from endocrine cells, inflammatory mediators from immune cells, 
neurotransmitter release from nerve cells, or the removal of waste molecules and bi- 
products. The /3-cells located in the islets of Langerhans in the pancreas are strongly 
secretory active in producing and releasing insulin. It seems that a better under- 
standing of these two processes could, e.g., support the early diagnosis of metabolic 
diseases, including diabetes mellitus type 2 (exocytic dysfunction, see Rorsman and 
Renstrom 1 8 9 1) or a more efficient delivery of insulin analogues (inducing endocy- 
tosis). 

Recent advances in observational and manipulative nano techniques and in meso- 
scopic coarse-grained computer simulation have provided substantial progress in 
visualizing, understanding - and possibly - influencing the bilayer membrane vesi- 
cle fusion (for a recent systematic review we refer to Shillcock and Lipowski 1 10 1). 
With present technology, however, the observational findings are hampered: geo- 
metric shape and the parameters of change can not be measured to the wanted de- 
gree of precision, simultaneously. So ffTOl, p. SI 196 deplores that "the molecular 
rearrangements that take place during the final stage of the fusion process, where 
the two initially distinct membranes join and produce a fusion pore, cannot yet be 
resolved by these experimental techniques" while "understanding how the stability 
of lipid membranes is overcome by the cellular protein machinery when required 
is a major topic of research" (I.e.). This is the challenge which we wish to address 
by summarizing various mathematical, electrodynamical and empirical facts. Seem- 
ingly, these facts have been overlooked or discarded in the literature. We shall show, 
however, to what extent these facts yield, in our perception essential, ingredients to 
a comprehensive understanding of this short very first phase of regulated exocytosis. 



2.2 Competing Mathematical Approaches to Space-Time Processes 

We shall distinguish three different mathematical-numerical approaches to model- 
ing spatial-temporal process of regulated exocytosis: highly aggregated compart- 
ment models, spatially distributed dynamical systems, and space-time integrating 
partial differential equations, where our focus will be. Clearly, all three approaches 
admit extensions from rigid, stiff and hence fragile deterministic to more robust 
stochastic modeling. Here, however, we shall not discuss such extensions. 

A first class of mathematical exocytosis models are Compartment Models, first 
introduced by Grodsky itTTI in 1972, assuming that there are two compartments 
(pools) of insulin granules, docked granules ready for secretion and reserve gran- 
ules. By assuming suitable flow rates for outflow from the docked pool and resupply 
from the reserve pool to the docked pool, the established biphasic secretion process 
of healthy j5 -cells could be modeled qualitatively correct. By extending the number 
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of pools from two to an array of six and properly calibrating all flow rates, Chen, 
Wang and Sherman in [ 12 1 obtained a striking quantitative coincidence with the ob- 
served biphasic process, see also Toffolo, Pedersen, and Cobelli [4| in this volume. 
The nice thing by such compartment models is that they invite the experimentalists 
(both in imaging and in proteomics) to verify the distinction of all the hypothetical 
compartments in cell reality and to assign biophysical values to the until now only 
tuned flow rates. A self-imposed limitation is the low resolution of the aggregated 
compartments which does not allow to investigate the local geometry and the energy 
balance of the secretion process. 

On the opposite length and time scale, we have a second class of mathematical 
models, namely the impressive numerical analysis of the bilayer membrane vesicle 
fusion by Molecular Dynamics (MD), Monte Carlo simulations (MC) and Dissipa- 
tive Particle Dynamics (DPD) on nanometer distances and fractions of nanoseconds, 
based on gravitational and electric forces between the particles, see the afore cited 
iflOl and Shillcock, [ 3 1 in this volume. Unfortunately, these computer simulations 
are also seriously hampered, namely by limitations of present hardware and soft- 
ware when one is addressing mesoscopic behavior, i.e., changes across many scales 
of the molecular characteristics - in spite of the impressive results when applying 
these methods to phenomena on the nano scale, like modeling the island dynamics 
of film crystallization in epitaxial growth driven by molecular beam epitaxy, see, 
e.g., Cafiisch and Li |fl~3). These limitations in present computer capacity require 
MD, MC and DPD simulations to make a priori assumptions about the pathway of 
the fusion process, e.g., spherical symmetry of the vesicles and planarity and circu- 
larity of the fusion pores - besides the often deplored "enormous gap between the 
sophistication of the models and the success of the numerical approaches used in 
practice and, on the other hand, the state of the art of their rigorous understanding" 
(Le Bris lfl4l in his 2006 report to the International Congress of Mathematicians). 
To keep these models transparent, a self-imposed limitation is their focus on the 
local neighborhood of the fusion event, neglecting long-distance phenomena like 
electromagnetic waves across the cell. 

Whereas compartment models respectively MD and DPD are built upon small 
respectively huge systems of ordinary differential equations with each unknown 
specifying temporal changes in one given pool or spatial box, we advocate a third 
class of mathematical models, namely modeling the dynamics and the geometry by 
partial differential equations. Consequently, we shall try to model the relevant pro- 
cesses by one or two spatial-temporal equations based on First Principles, instead of 
the few aggregated purely temporal pool equations in compartment modeling or the 
three millions of purely temporal equations for spatially distributed boxes in ifTUll . 
p. 1 197 (which still is a very poor particle number for a 100 nm x 100 nm x 42 nm 
simulation box). 

Moreover, the simplicity of our fundamental equations admits a transparent in- 
corporation of long-distance phenomena. In such a way, our approach takes its point 
of departure not only in the rather well studied elastic and electric properties and 
potentials of and across the plasma membrane and the viscosity of the cytosol, but 
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in non-stationary, dynamic electromagnetic properties. To us, the basic electromag- 
netic character of the fusion process becomes evident in 



• the observed electromagnetic (wandering) field waves, see Sect. 2.3 

• the closure of the corresponding magnetic wave over the plasma membrane, see 
Sect. [24] 

• the observed forming of a narrow dimple, solely to be understood like a capaci- 
tator, see Sect. |2.5[ and 



the observed flickering of the secretion process corresponding to the natural vari- 



ability of the (wandering) field wave generation, see Sect. 2.6 
The goal of our approach is 

1. to develop a simple free boundary model for the dimple forming process, see 
Sect. [3] 

2. to focus on regularity and possible singularity of the free boundary, see Sect. [4] 

3. to provide a reliable framework for estimating (and, hopefully, influencing) the 
parameters which control the speed of the process, see Sect. |5.2| 

4. and to formulate a bundle of model-based observation plans to verify or falsify 



our assumptions, see Sect. 5.3 



Our approach is inspired by recent work of A. Friedman and collaborators about 
tumor growth, see lfT5l - lfl8l : by a theoretical analysis of adhering lipid vesicles with 
free edges, see [ 19 1 by Ni, Shi, and Yin; by the electrodynamic challenge to under- 
stand the observed Ca 2+ oscillations, rightly perceived as being "contradictory and 
often do not support the existing (electrostatic) models" (Fridlyand and Philipson 
[20]) and by the mathematical challenge to understand the (in elastic terms counter- 
intuitive) dimple formation so well described in the literature, see, e.g., Monck and 
Fernandez [23], Rosenheck l24l . Lentz et al. [7|, Koch et al. 0. 



2.3 Oscillatory Intracellular Release and Binding of Ca 2+ Ions 

We recall a few basic observations of Ca 2+ -oscillations and postulate a simple but 
powerful procedure of generating an electrodynamic field. 



2.3.1 Basic Observations of Ca 2+ -oscillations 



By fluorescence microscopy, empirical evidence has been provided about pulsating 
Ca 2+ activity at extreme low frequency / ~ 0.1 Hz <C 3 Hz (for comparison, the 
house low frequency grid is of 50Hz, i.e., spikes in intervals of 20 msec), prior to the 
fusion event, see Kraus, Wolf and Wolf, ll25ll and Bernhard Wolf's homepage [26| 
with informative video animations of calcium oscillations, and the comprehensive 
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review and analysis by Salazar, Politi and Hofer, ll27ll {j Wolf's observations were 
made with HELA cancer cells. Hofer et al. gives a general model, based mainly on 
observations in muscle cells. 

For corresponding observations for j3 -cells, we refer to Maechler [28 1 and 
Fridlyand and Philipson [20| in this volume, who deal with various types of low- 
frequent oscillations. 

The following can be seen in many cell types: when a cell is stimulated ("pol- 
luted") by a Ca 2+ -mobilizing stimulus, the changes in the cytosolic calcium con- 
centration [Ca 2+ ] c occur as repetitive spikes that increase their frequency with the 
strength of the stimulus (see also Berridge, Bootman, and Roderick ll30l and Gaspers 
and Thomas OTl ). It is well known that an increase in [Ca 2+ ] c , ultimately, regulates 
a plethora of cellular processes mediated by Ca 2+ -dependent enzymes that, in turn, 
modify downstream targets commonly by phosphorylation. 

Investigating Ca 2+ decoding in an analytically tractable generally applicable 
model, Hofer and collaborators address the question "Under which conditions are 
Ca 2+ oscillations more potent than a constant signal in activating a target protein?" 
in 11271 p. 1204] in complex biochemical terms of binding and release rates of Ca 2+ 
ions. 



2.3.2 Postulated Electrodynamic Field Character 

We wish to supplement these investigations by recalling pieces of circumstantial 
evidence which may support the hypothetical electromagnetic character of the Ca 2+ 
oscillations. As said before, Ca 2+ handling is extremely complex. 

1. There are various Ca 2+ -storage organelles in the j3 -cells, first of all the spatially 
distributed and clearly separated mitochondria and the spatially rather extended 
smooth endoplasmic reticulum - SER, see also ll20l Section 5]. 
item There is some agreement in the literature that the handling by the mi- 
tochondria, and not the SER, is decisive for, at least, the slow Ca 2+ oscilla- 
tions: "Metabolic profiling of /3-cell function identified mitochondria as sensors 
and generators of metabolic signals controlling insulin secretion", according to 
Maechler |28l : and Fridlyand and Philipson EDI refer to "data demonstrating 
that slow oscillations can persist in the presence of thapsigargin, the agent that 
blocks SERCA and empties the ER stores...". 



1 Our model cell is a pancreatic j3-cell where a single release is slow and may take seconds. Cor- 
respondingly, we expect a low Ca 2+ oscillation rate in intervals in the range of seconds yielding 
extreme low frequency of the observed 0. 1 Hz. For nerve cells the reaction time, and so the release 
time, is in the range of msec, possibly less than 100 (isec, see Jahn, Lang, and Siidhoff [29 |. Corre- 
spondingly, we expect a high Ca 2+ firing rate in intervals of, e.g., 10 msec yielding a frequency of 
100 Hz with associated high energy losses. So, our electromagnetic free boundary route to vesicle 
fusion can not function in nerve cells unless the neurotransmitter vesicles are kept waiting very 
close to the plasma membrane. Here our finding coincides with the well-known deviating high 
energy consumption of nerve cells. We shall elaborate this aspect in Sect.[5]below. 
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2. To us, these observations fit nicely with our postulate of the electrodynamic char- 
acter of the slow Ca 2+ bursts: spatially distributed and temporarily coordinated 
(not necessarily simultaneous) release and uptake of Ca 2+ ions can generate an 
electrodynamic field with magnetic character. 

Of course, there are many mechanisms and systems that contribute to the release 
and uptake of Ca 2+ from intracellular organelles, some of these even have bell- 
shaped effects on Ca 2+ release, and phenomena such as Ca 2+ -induced Ca 2+ release 
also take place (these facts were communicated to the authors by Pociot and St0rling 
||2T)). In the following, we anticipate the existence of a (not yet fully confirmed) 
system of aggregated build-up of an electrodynamic field with magnetic character 
in /3-cells upon glucose stimulus: In short, we have movements of Ca 2+ ions in and 
out of the mitochondria. Movements of ions are currents. In-and-out movements 
are AC currents. AC currents produce AC fields, or rather, AC field oscillations. 
Superposed oscillations may produce moving fields. 

In this way, we suppose, the organelles build an alternating electric current den- 
sity (also "displacement vector field") D of low frequency by superposition of spa- 
tially distributed, temporally coordinated and directed Ca 2+ activity. As usual in 
electrodynamics, we shall speak of two different electric fields, D and E. This sec- 
ond electric field E is given by the relation D = eE, where e = EoE r denotes the 
dielectric constant. Note that our writing of all electromagnetic units and equations 
follows Jackson ll22l in the units V, A, sec, and m with, e.g., kg=VAsec 3 /m 2 of the 
System International - SI, which is predominant in engineering literature. 

A physical model of generating electromagnetic field waves by spatially and tem- 
porally distributed excitation was built by Koch and Stetter, see |http: //www.| 
|f eldkraf t . de7| It is called Dynamical Marker and consists of a couple of coil 
arrays, electronically regulated for direction-, frequency- and amplitude-controlled 
generation of a field wave. The instrument has been applied in various cell physi- 
ological experiments for slow and efficient transport of beads into cells across the 
plasma membrane. 

There are a couple of related questions which will require separate investigation 
(see also below Sect. |5.3| >: 

1. How strong is the evidence that the Ca 2+ oscillations, contrary to the membrane 
process of Ca 2+ bursts, originate from an array of Ca 2+ depots (SER and mito- 
chondria), organized in directed lines? 

2. How does the cell select the Ca 2+ storages to participate in the generation of the 
alternating current? 

3. How is the sequential release and binding of the Ca 2+ ions of the different stor- 
ages controlled, i.e., how is the correct spatial and temporal coordination of re- 
lease and binding obtained? 

4. What role plays the observed branching of mitochondria in active j3 -cells, con- 
trary to the dipole shape of mitochondria in tired j5 -cells? 

5. Can the magnetic character of the field wave, which is produced by low frequent 
Ca 2+ oscillations, be influenced via an external field with similar character? 



Chapter 30 



- Geometric and Electromagnetic Aspects of Fusion Pore Making 



B 



11 




Fig. 2 Functional heterogeneity: left, a multitude of branched mitochondria in a vigorously re- 
sponding /3-cell; right, relatively more non-branched dipole mitochondria in a less active cell at a 
comparable locus; branch points are highlighted with red spheres. The panels are from Noske et 
al-> (32 1 and reproduced with permission of Elsevier. High resolution originals were supplied by 
courtesy of B. Marsh, University of Queensland, Brisbane, Australia. 



To Questions 1 and 2, it may be mandatory to distinguish between two different 
types of Ca 2+ burst: the Ca 2+ oscillations addressed here are prior to the secretion 
and are, presumably, generated by arrays of activated calcium depots distributed 
through the full length of the cell. On the contrary, the Ca 2+ influx through the ion 
channels of the plasma membrane during regulated exocytosis is mainly effective 
close to the plasma membrane where it increases the concentration [Ca 2+ ] t and 
changes the electric potential across the plasma membrane. It is worth mentioning 
that the array character is evident from the observed oscillations in the form of 
directed, oriented waves. 

To Question 3, we imagine that the coordination of the activity of the partici- 
pating Ca 2+ storages is not controlled externally, e.g., by the nucleus, but happens 
spontaneously by selforganization: we notice that the storages of molecular calcium 
sense and respond to stimuli by periodic release and binding of Ca 2+ ions and sup- 
pose that the sequential coordination between spatially distributed loci of release 
and binding minimizes the energy consumption for maintaining the activity and for 
establishing a suitable average concentration [Ca 2+ ] t . 

To Question 4, we recall from Noske et al [32] a remarkable discovery. The 
authors imaged and reconstructed two j3 -cells from the same glucose-stimulated 
mouse islet by single axis, serial section electron microscope tomography (ET) at 
magnifications of 4700 x and 3900 x, respectively, that resulted in whole cell to- 
mograms with a final resolution of 15-20 nm. In addition, they developed several 
new methods for the abbreviated segmentation of both cells' full complement of 
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mitochondria (i.e., the most prominent Ca storages) and insulin secretory granules 
for comparative analysis. 3D reconstruction by ET of each of the two j3 -cells (des- 
ignated 'ribbonOl' and 'ribbon02') indicated that ribbonOl responded more vig- 
orously to glucose-stimulation than ribbon02 and contained about twice as many 
branched mitochondria (26 out of a total number of 249 mitochondria) as ribbon02 
(with 10 branched mitochondria out of a total of 168 mitochondria). See also the 
recent Marsh and Noske [33 1 in this volume. 

Now, from our electromagnetic point of view, the advantage of branched mito- 
chondria is clear for the generating of a field wave (i.e., the pulsating Ca 2+ oscil- 
lations through the whole length of the cell): a single (non rotating) dipole can not 
generate or initiate a (directed, oriented) field wave. That requires a branched struc- 
ture with spatially and temporally shifted serial activity, as demonstrated also in the 
design of the mentioned simple field generator by Koch and Stetter. Of course, na- 
ture's regulation of the ion firing may be much more sophisticated than the crude 
engineering design of the "Dynamic Marker". After all, the eukaryotic cells had 
many more years to test and optimize different designs in evolution. 

To Question 5, we refer to an experimental setting described below in Sect. |5.3| 



2.4 The Magnetic Character of the Induced Field Wave 

Perhaps one of the most delicate of Maxwell's equations (cf. Box) is his modifica- 
tion of Ampere's Law by adding the displacement current density on the right, 
i.e., the electric side: 

curl H = J+^, respectively J Hds = J J (J+^jdA. (1) 

Here H denotes the oriented, electrically generated magnetic field, B = /Xo/X,H the 
corresponding magnetic flux density (also "magnetic induction") of frequency / 
and amplitude B, and J the current density vector for a conductor crossing an area 
A which is bounded by a contour C. Moreover and jj,, denote the absolute and 
relative magnetic permeability. 

Equation ([T} is our basic equation regarding the magnetic character of the ob- 
served Ca 2+ oscillations. Note that our field wave can not be compared with an 
electromagnetic high frequency wave in radio transmission. Its propagation veloc- 
ity is comparable to sound waves in water and far below light velocity. Moreover, 
because of the low frequency of our oscillations, the displacement current density 
is relatively small. That indicates that the magnetic character dominates the elec- 
tric character of the field wave. 

While a wave with dominant electric character has large losses in cytosol (which 
is comparable to salty water in its electric conductivity), the magnetic character 
makes the propagation almost free of losses. To sum up, by the described intra- 
cellular Ca 2+ oscillations a field is generated with only marginal losses because 
the transmission is almost independent of the material constant £,. . As a result, the 
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The four Maxwell's equations at a glance 

Gauss' law V ■ D = p. describes how an electric field is generated by electric charges. 

Gauss' law for magnetism V ■ B = valid for the absence of magnetic charges (also called 
magnetic monopoles), analogous to electric charges. 
dB 

Faraday's law V x E = describes how a changing magnetic field can create (induce) an 

dt 

electric field. 

Ampere's law with Maxwell's correction VxH = J , + states that magnetic fields can be 

' dt 

generated in two ways: by electrical current (this was the original Ampere's law) and by changing 
electric fields (this was Maxwell's correction). 

E is the electric force field, measured in Volt per meter, and 

D is the electric flux field or displacement field, also called the electric induction or the electric 
flux density. It is measured in Coulombs per square meter. 

H is the magnetic force field or magnetizing field, also called the auxiliary magnetic field, the 
magnetic field intensity or just the magnetic field. It is measured in Amperes per meter, and 
B is the magnetic flux field, also called the magnetic induction, the magnetic field density or the 
magnetic flux density. It is measured in Tesla. 

p r is the free charge density, measured in Coulombs per cubic meter. 
J. is the free current density, measured in Amperes per square meter. 

More to read: J. D. Jackson, Classical Electrodynamics, 3rd ed., John Wiley, New York, 1999. 

Added by the editors 



moving Ca 2+ ions from the intracellular distributed storages are in fact an AC cur- 
rent, generating a nearly loss-free moving (wandering) magnetic field wave which 
transfers energy to a selected transmembrane subregion of the cytosol between the 
plasma membrane and a single vesicle. As we shall explain, this energy and the 
corresponding forces act on the free ions distributed in one or other way among the 
phospholipids of the plasma membrane and pull the plasma membrane towards the 
vesicle. See also [34] for a video animation of a macroscopic field wave. 

Our biophysical approach is classical, hence we assume that there are no sources 
for magnetism, no magnetic monopoles, at least not present in our /3-cell. Con- 
sequently, we obtain from Maxwell's equations divB = 0, i.e., the magnetic field 
wave induced by the alternating current is closed in the sense of vector analysis and, 
consequently, the path of the field wave is closed. 

Of course, it must be investigated in detail how the magnetic wave is closed. 
Fluorescence microscopy gives the impression that the observed Ca 2+ oscillations 
are collective phenomena of cell ensembles: the oscillations propagate through the 
ensemble like chained waves. We, however, assume that all magnetic waves are 
separated from each other and are closed over the single plasma membranes. One 
reason is that the plasma membrane is perforated by a multitude of ion channels ere- 
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ated and maintained by the presence of enzymes like various kinases and phospho- 
lipases. Most encymes contain Fe atoms, see Jensen 1331 pp. 134f] and the review 
on iron biominerals |36|. Consequently, the magnetic field wave will search for a 
circuit through the plasma membrane. 

While we have put the determination of the Fe content of the plasma membrane 
on our experimental agenda below in Sect. |5.3| we should mention that the mag- 
netic permeability ji r can not be measured directly. However, there are methods to 
determine the magnetic permeability ji, ~ 1 .0000007 of hemoglobin of deoxidized 
venous blood non invasively in a precise way. Similar methods will be applicable 
for investigating the plasma membrane of life cells. Roughly speaking, the indirect 
methods work by comparative measurements after inflicting a magnetic pollution 
on a harmonic oscillator. 

Now it is not difficult to understand the making of the fusion pore and the dim- 



ple formation (see the following Sect. 2.5 for the empirical evidence) in qualitative 



terms. Let us fix the notation. By electron microscopy, we can distinguish the fol- 
lowing clearly separated regions: 

Do Amorphous outside cell neighborhood 

D\ Plasma membrane with boundary dD\ = r]_o U F[_2 

D 2 Cytosol 

Do, Vesicle membrane with boundary dD?, = i~3_2 U I3-4 

D4 Vesicle lumen 

Mj Activated molecular Ca storage organelles, j = 1 , . . . n (not depicted) 

N Cell nucleus (not depicted). 



Fig. 3 Plasma membrane 
section with cone like dim- 
ple (left) under transversal 
displacement u and spher- 
ical vesicle (right) before 
docking. For a regular cone 
shaped dimple, the two 
bold dots mark the circular 
base line r(t) = d{(x,y,z) \ 
u(x,y,z,t) > 0} of the dimple 
at time t. In space-time, the 
union {T(f) x {t}},>o of the 
t -components F(f) form the 
free boundary F, see below 
before Eq. \\6) and 



Sect. 
Sect. 



4.1 



4.4 
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Note that the region D\ (the plasma membrane) consists of a phospholipid bi- 
layer. It has a surface il_o as its outside boundary (towards the amorphous outside 
cell neighborhood Do) and a surface Ij_2 as its inside boundary (towards region 
£>2 of the watery cytosol within the cell close to the plasma membrane). Moreover, 
we have the region D3 consisting of the vesicle membrane (same material like Di) 
and the region D4, the interior of the vesicle, containing the material to be released 
through the plasma membrane. Moreover, there is a multitude of activatable Ca 2+ 
storage organelles {Mj} spread through the interior of the cell. Finally, we have the 
cell nucleus N, see the abstraction of Fig.[T]in Fig. [3] Typical approximate diameters 
are 10 pim for most animal cells; 5 fim for the nucleus; 100 nm for the vesicles; and 
8 nm for plasma and vesicle membrane. 

In the preparation of the fusion event and the making of the fusion pore, there are 
apparently only two active regions, in addition to the cell nucleus and the mitochon- 
dria and other Ca storages, namely the plasma membrane D\ and the cytosol Dj_. 
The plasma membrane D\ forms a conical inside oriented dimple (pit) towards the 
vesicle of around 10 nm base diameter and 10-20 nm height. For the true lipid bi- 
layer membrane-vesicle fusion event, following the making of the dimple, so-called 
transmembrane proteins become active in the cytosol region D2 between plasma 
membrane dimple and the vesicle and pull and dock the vesicle membrane D3 to the 
plasma membrane D\ over a distance of up to 100 nm (see also Fig . [TJ . 

When the deformation of the plasma membrane is sufficiently sharp it ends in 
a branch point, i.e., a singularity of Ij-2> called hemifusion. Then it comes to a 
break-through (called fusion pore), and the content of the vesicle begins to diffuse 
from the vesicle compartment D4 into the outside region Do. It appears that this 
process sometimes is interrupted (so called flickering, see below), i.e., the fusion 
pore is hardly maintained by elastic forces alone but needs probably the presence of 
an electromagnetic field and is interrupted when this magnetic field is interrupted. 

What is controlling the well-functioning of the fusion event? 

Working hypothesis 1: The regions 

D (X c = 0), Di (variable X c (x,t)) and D 2 (X c = 0) (2) 

are distinguished by their capacitive reactance Xc := 1 /((oC) where co — 2nf with 
constant / and C denotes the capacitance. Note that forming the dimple produces an 
increase of the dielectricum (between the "plates") and so implies increasing C and 
decreasing Xc until Xc vanishes in the fusion pore. 

Working hypothesis 2: The vesicle is densely packed with material and so not 
subjected to deformations easily. 

Working hypothesis 3: We envisage the following feedback mechanism for 
forming the dimple and preparing the fusion event. 

1. The Ca 2+ ions from locally distributed intracellular Ca-molecule storages start 
low frequent oscillations, as described in 11251 . which are superposed in a con- 
trolled way, and a dynamic field wave is produced pointing to a specific region 
D2 jCr it selected for most suitable membrane-vesicle fusion. In the beginning, the 
magnetic flux density vector B ("magnetic induction") is low because the mag- 
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netic wave does not easily enter the plasma membrane D\ to close itself in a 
circuit because of the originally high Xq in D\ . 
2. The form change decreases X c close to the emerging dimple. That permits the 
magnetic wandering wave to enter D\ more easily and so increases its current 
density (the sharpness of its pointing) and its amplitude B. And so on. 




Fig. 4 Cross section through plasma membrane with dimple and, at the bottom, a glimpse of the 
granule. Electron micrograph supplied by courtesy of D. E. Chandler, Arizona State University, 
Tempe, Arizona. 
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In the introduction to this chapter we defined the fusion pore as the molecular struc- 
ture that transiently connects the lumens of two membrane compartments during 
their fusion. We emphasized that making the fusion pore plays a key role in all 
intracellular trafficking and endocytotic and exocytotic pathways in all eukaryotic 
cells, including the regulated exocytosis in endocrine, exocrine and neuronal cells 
like our j3-cell. However, from Monck and Fernandez 11231 7992] to Shillcock and 
Lipowski [10, 2006], researchers agreed that despite its importance, the nature of 
the fusion pore is unknown ( 1231 p. 1395]). 

In a remarkable series of micrographs, based on rapid freezing techniques for 
electron microscopy, the renown expert on mammal egg cells, Douglas Chandler 
from Arizona State University and collaborators captured the formation of the fu- 
sion pore in mast cells already some 30 years ago, [37, 38 , 39]. They demonstrated 
that the pores are made of a curved bilayer which spans the granule and plasma 
membrane. The micrographs also gave a hint of the events preceding the making 
of the fusion pore: namely the formation of a dimple that approaches the granule 
membrane after stimulation, see Figures [4] and [5] 

Since then, the dimple formation has been observed in many different cell types 
immediately upon stimulation before the fusion event, see, e.g., Jena and collabo- 
rators [40 1 . Using atomic force microscopy (AFM) they demonstrated the presence 
of many simultaneous dimples in pancreatic acinar cells after exposing them to a 
secretagogue. The paper contains references to analogous demonstrations revealing 
the presence of pits and depressions also in pituitary and chromaffin cells prior to 
secretion. To the best of our knowledge, capturing the very fusion event of regu- 
lated exocytosis in pancreatic j3 -cells has not yet been achieved by present imaging 
methods. 

Not surprisingly, the secretory granules have never been seen to form dimples 
on their own membrane, in accordance with the common perception that the plasma 
membranes are relatively slack and the membranes of densely packed granule under 
tension. That corresponds to our preceding Working Hypothesis No. 2. 



2.6 The Flickering of Regulated Exocytosis 

Another feature of the fusion pore making requires explanation (and is nicely ex- 
plained by instabilities of the discussed AC current), namely the flickering of the 
fusion event, i.e., the common observation that the fusion pore can be maintained 
only after a while of opening and is not stable immediately after its making. Irregular 
rapid pore openings and closures are observed that last from a few milliseconds to 
many seconds, see Fernandez et al. ETll for fusion-pore flickering (kiss-and-run) in 
mast cells, and Rosenheck [ 24 1 and Jahn et al. 1 29 1 for wider reports on the observed 
flickering of the fusion event, mostly for synaptic vesicle exocytosis. 
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Fig. 5 The fusion pore shown as imprint of the plasma membrane dimple onto the granule mem- 
brane, from 1 39 Fig. 9C]. Reprinted with permission by Springer- Verlag. The high resolution mi- 
crograph was supplied by courtesy of D. E. Chandler, Arizona State University, Tempe, Arizona. 



3 The Model 

In this section we fix our notation and introduce the basic equations for the propa- 
gation of the electromagnetic (wandering) field wave and the making of the dimple. 



3.1 The Force Balance Equation 



Let r = r(x,y,z,t) denote the displacement vector of the dimple and m the mass of 
an elementary unit of the dimple. Then the resulting force for making the dimple is 

approximately equal to m^-ir, i.e, we can write 

d 2 r 



m ~fo2 = F elas+F vi i- 



(3) 
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where F e i as denotes the restoring (or elastic) force, F V7S stands for viscosity reac- 
tion of the medium surrounded our plasma membrane, and F ext is an external force 
providing the membrane displacement from equlibrium state. 

Our first task is to describe the restoring force F e ] as . In first approximation, we 
shall assume that the plasma membrane is a surface in M? without bending resis- 
tance. Hence F e i as is defined only by the surface tension and the variation of the 
membrane surface area. 

To describe F e [ as and F v ; s more precisely, we adapt the standard text book 
model for the suspended vibrating string, respectively, vibrating plate (see, e.g., 
Churchill [42, Sect. 93] and Logan l43l ). We consider ID- and 2Z)-membranes 
separately, because the restoring force behaves differently in one-dimensional and 
many-dimensional cases. We begin with the ID-case, since all the arguments are 
simpler in that situation. 



In the equilibrium state, let our ID-membrane coincide with the x-axis; let u = 
u(x,t) denote the displacement of our plasma membrane from equilibrium state at 
the point x and at the time t ; and let p = p (x) denote the linear membrane density at 
the point x. We restrict ourselves to sufficiently small deformations, so, for now, we 

du 

will neglect all the terms that are of higher infinitesimal order with respect to 



3.2 The ID-Case 



u ZA 




T(x+Ax) 



x 



a(x+ax) 



Fig. 6 Tension along the plasma membrane 
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Since our plasma membrane has no bending resistance, its tension T(x,f) at the 
point x and at the time t is directed along the tangent to the membrane at x. There- 
fore, the unit (x,x + Ax) is subjected to tensions T(x + Ax,t) and —T(x,t) (see Fig. 
[6j. Moreover, according to Hooke's law, |T(jc,f)| does not depend on x and t, i.e., 
|T(*,f)|=T . 

Defining T\ as a relaxation time due to the action of the surrounded medium, we 
note that F v i s is directed parallel to the vertical coordinate axis and can be modeled 
1 du 

as being proportional to =— . 

T\ at 

Let us denote by f(x,t) the density of an external force F ext acting on the mem- 
brane point x at the time t and directing along the vertical axis. For a description of 



our electromagnetic candidate, see below Sect. 3.5 

Thus, projecting Eq. <(3j onto the vertical coordinate axis and taking the preceding 
relations into account, we get the equality 



pAx 



T(x + Axj)sm(a(x + Ax)) - T(x.t)sma(x) - —Ax^" 1 +f(x,t)Ax. (4) 

7] at 



In the context of our approximation 

tan (a) 



du 



sin(a) = — = » tan (a)- , 

v /l+tan 2 (a) ox 



and, consequently, equality (Hh takes the form 



d 2 u To 
dt 2 pAx 



du(x + Ax,t) du(x,t) 
dx dx 



2 du 1 , s 



Passing in <|5j to the limit as Ax — > we arrive at 



d 2 u 2 du 2 d 2 u 1 
-dfi + T l J t - C >dx 2 - = p f{x > t) 



To 



where cf = — is the speed of sound in the membrane. 
P 



(5) 



(6) 



3.3 The ID-Case 

Similarly to the ID-case we assume that in equilibrium state our plasma membrane 
lies in a subspaceXF (see Fig. |7j> and u = u(x,y,t) denotes the membrane displace- 
ment from equilibrium state at the point (x,y) and at the time t. We will consider 
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only small deformations such that 



du\ 2 ( <9m x 2 



dxj <<h [Ty) <<L (?) 

Let da be the unit length of some closed path lying on the membrane surface, 
and let P be a point belonging to da. Then the unit da is subjected to the tension 
force Ida, where T = T(x,y,t) denotes the surface tension. Due to absence of the 
membrane resistance to bending and shear, we can say that the vector T always lies 
on the hyperplane L\ tangential to the membrane surface at the point P, and T is 
orthogonal to da (see Fig. [Til. 

In addition, inequalities (|7]i guarantee that the tangential hyperplane L\ lies al- 
most parallel to the hyperplane XY. To prove this statement, we have to show that 
the length of the projection of the vector T(x,y,t) onto XY is approximately equal 
to |T(x,y,f)|. Indeed, by definition 



|Projection XF T(x,y,f)| = |T(x,y,f)|cos(j3), 

where j3 stands for the angle between the tension vector T and hyperplane XY. It is 
easy to see that j3 is not bigger than the angle 7 between the tangent hyperplane L\ 
and XY. Therefore, 
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COS Q3) ^ COS (7) = 
and, consequently, 



-1/2 



1. 



|Projection zy T(x,y,f)| w |T(jc,y,f)|. 

Hooke's law guarantees that |T(x,y,?)| does not depend on the f-variable, whereas 
the orthogonality of T(x,y,t) and da provides the independence of |T| on the vari- 
ables x and y as well. It means that 

|T(jc,y,?)| = To = const. 

Now, considering a rectangle unit ds = AxAy on the membrane surface, we can 
write the restoring force acting on this unit as: 



T Ay 



■■T Ay-^Ax- 



-T Q Ax 



du 
dy 



du 



d 2 u 
ToAx^Ay: 



d*u 
dx 2 



d^u 
dy 2 



AxAy. 



It remains to describe the external and viscosity forces acting on ds. Similarly 
to the ID-case, f(x,y,t) denotes the density of the external force F ex t at the point 
x and at the time t. It is directed orthogonally to the membrane surface, while F v i s , 

du 

directed opposite to the vector F ext , is proportional to the velocity . 

dt 

Let p(x,y) denote the membrane surface density, then the mass of the unit ds 
is equal p(x,y)AxAy. Finally, defining a relaxation time due to the action of the 
surrounded medium by T\ , we can write the variant of Eq. Q for 2D-membranes as 
follows: 



d 2 u 

P AxAy^=T 



d^u 
dx 2 



d^u 

dy 



2p du , . 

T\ AxAy ~di +f( x >y> t ) AxA y- 



After elementary transformations, Eq. ((HI) takes the form 



dhi 
dt 2 



2 du 



d 2 u d 2 u 
dx 2 dy 2 



-f(x,y,t), 

r 



c; 



(8) 



(9) 



3.4 Further Approximations 



Having equations (|6) and |9) at hand, we observe that the process of dimple forming 
is quasi-static, i.e., 
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d 2 u 

The latter means that we can neglect this term in both equations. 

It should also be pointed out that u = \Du\ = at those membrane points where 
there is no influence of the external force. Here Du denotes the spatial gradient of 
the displacement u. 

Recall that the characteristic function %e of a set E is defined by 



Xe(z) = 



1, forze£\ 
0, forz^£. 



Taking into account all the above remarks we get the following model equations 
for the dimple forming in the one- and two-dimensional cases, respectively: 



2 du o d u 1 



7i dt ° s dx 2 p 
2 du 2 ( d 2u d 2 u\ 1 
Yi ~dt ~' 



= -f(x,t)X{ u >o}, (10) 



(d 2 u d 2 u\ 1 . 
W + dy 2 -) = p fiX ^ t)X{u> ° } - 



3.5 Lorentz Force 

Clearly, there is a variety of external forces resulting from the electromagnetic field 
wave. One can expect that they all play together in forming the dimple. However, 
taking our point of departure in an alternating current, we shall concentrate on the 
Lorentz force, (implicitly deal with the Coulomb force), but discard the dipole elec- 
tric force, the magnetic force, and the van der Waals force for now. 



3.5.1 Peculiarity of the Lorentz Force 

In our physiological context, the peculiar role of the Lorentz force is that it exerts its 
action in one fixed direction, the direction of the propagation of the field wave, even 
when it is related to an alternating current. Roughly speaking, it is the sophistication 
and power of the electromagnetic aspects of the regulated exocytosis that a relatively 
weak and extremely low frequent electrodynamic wave can transport energy along 
a straight line over a large intracellular distance and exert its action on the charged 
phospholipid molecules in the plasma membrane. These charged molecules make 
a substantial part of the eukaryotic plasma membrane, around 11% according to 
Rosenheck [24 1. 
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High voltage devices: quantitative comparison of electric field strenghts in electrical power 
plants and animal cells 

A conventional 36.5 MegaWatt, 120 rotations per minute, 6.600 Volt High Tension Synchronous 
(HTS) generator yields a field strength of 100.000 V/m. 

For comparison, in animal cells we consider the electrostatic field over the plasma membrane of 7 
nanometers diameter and a potential of about 70 milliVolt. That yields a field strength 



\E\= 



70-KTV 



ao 7 



7 10" s m 

which is a 100 times stronger than the field strength generated in an electric high power station. 

Note that electric power engineering has to handle electromagnetic fields oscillating with low and 
very low frequencies. Therefore, in power engineering, the electromagnetism can be split into its 
electric and its magnetic character for optimizing the design of large machinery. Low frequent 
dynamic oscillations can yield (travelling) field waves, an aspect interesting for biophysics. 

More to read: M. Koch, http://www.freldkraft.de/ 



3.5.2 Energy Estimates 

It may be instructive to have a rough idea of the scale of the electric forces and 
related energies around the dimple formation, see also the box. Based on the esti- 
mates given by Rosenheck I.e., we have a surplus of about 10 10 (negative) charges 
per 1000 nm 2 plasma membrane area, i.e., about 10 9 charge carriers around the 
dimple top stretching over an area of approximately 100 nm 2 . The energy of 
one charged phospholipid molecule was calculated by Rosenheck I.e. as being 
10~ 19 Wsec = 10~ 19 Nm. Dividing by a characteristic length of lOnm we obtain 
a force of around 10~ n N, i.e., a total force of 10~ 2 N exerted on the dimple region. 

For comparison, the gravitational force of the mass 10~ 21 kg of a dimple of about 
10 3 nm volume (and specific weight comparable to water) would only sum up to 
10" 21 x 9.8 - 10- 20 N. For another comparison, we refer to the Koch-Stetter elec- 
tromagnetic field wave generator. For now, its B-field is of about 35 mT (milliTesla) 
and the exerted forces are visible and can be measured. 

To sum up, for the making of the dimple we are investigating electric forces 
which surpass by large the above estimated 10~ 2 N. 



3.5.3 Model Quantities 



Our first model quantity is the capacitive reactance Xq as scalar function of place and 
time, depending on the shape of the dimple. Our second model quantity is the time 
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dependent (namely dimple shape and Xc dependent) amplitude B of the magnetic 
flux density. 

3.5.4 Model Equations 

As mentioned before, we have two model levels. 

(I) There must be a variational equation, minimizing an energy functional or 
another related cost functional which gives the change and spatial distribution of 
X c - 

(II) As discussed above, there are also elastic and viscous forces resisting the 
necessary re-packing of the lipid heads under deformation and pulling the plasma 
membrane back in the more smooth non-dimple form, interrupting exocytosis when 
the electrical and magnetic production is interrupted. 

3.5.5 The Lorentz Force 

Our favorite external force for making the dimple is the Lorentz force 

F L = < ? E + < ? (vxB)-7v, (12) 

where q is the charge (i.e., twice the number) of the released Ca 2+ ions, Eo is the 
background electric field, here assumed to be zero, B is the magnetic flux density, 
v is the velocity of the charged dimple particles with respect to the coordinate sys- 
tem in which the quantities F, Eo, and B are calculated, i.e., v is up to the sign the 
wandering velocity of the directed electromagnetic field wave, and yv is the friction 
force under propagation. Once again, equation ([12} is written in SI. In the cgs sys- 
tem, common in physics, the term q{\ x B) has to be divided by c, the velocity of 
light in vacuum. 

Note that we have for a (moving) field wave (see, e.g., [44, Chapter 13]) 

B(x,t) =Bcos(k-x(f)-Jtv(x,f)f) +B DC , (13) 
E(v)=v(v)xB(v). (14) 

Here, Bdc denotes the background B-field corresponding to the direct current Eo , 
x = (x,y,z) denotes a position, k denotes the wave vector with k = |k| and v(x,f) = 
|v(V)|. 

The Lorentz force of Eq. ( fl2] > can be inserted into the general balance equation 
([3]), into the 1 -dimensional model equations Q, ([6}, simplified to (jTUJ, and into the 
2-dimensional model equations (|8), simplified to ( flT) . 
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3.5.6 Work Equation 



In Sect. 2.4 we have sketched the feedback mechanism of the dimple forming and 
explained why and how the forming of the dimple strengthens the electromagnetic 
field wave. It is beyond the range of this chapter to model that mechanism in detail as 
a free boundary problem. Basically, to relate our defining quantity Xq to the listed 
balance equations, we shall express the power in terms of Xc and then derive an 
integral for the electromagnetic energy density (per volume), where Xc enters. The 
details will be worked out separately. 



4 Apposite Results on Parabolic Obstacle Problems 

At the end of Sect. 3.4 we introduced the model equations ( fT0| ) and ( fTT) that can be 



treated as parabolic free boundary problems (FBP). 



4.1 Review of Free Boundary Problems 

The expression free boundary problem means that we deal with a problem with 
two a-priori unknown objects: an unknown set coming up in a solution of a partial 
differential equation. 

A typical example is the Stefan problem describing the melting of an ice cube in 
a glass of water. If ice begins to melt then the region occupied by water will grow 
and the interface-surface between the ice and the water (it is called the free bound- 
ary) will move and change its shape, see, e.g., Friedman B31 Sect. 1.9]. As another 
typical example of FBP we mention the flame propagation problem describing the 
evolution of the flame front, for that see the various contributions in [46 , Chapter 8]. 

Using a common transformation of the independent variables, we may normalize 
the coefficients and so reduce our model equations (jTOf and ( fTTj ) to the following 
problem: 

Au(x,t)-~(x,t)=f(x,t) X{u>0] , almost everywhere (a.e.) in ^, (15) 
u(x,t) ^ 

where (x,t) denotes the points in R" x R with the space variable x = (jci , ,x n ) 

belonging to R" and the time variable t belonging to K, A is the Laplace operator 
defined as 

d 2 u 



AU = ^dx^ 

1=1 OX i 
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X{u>o} is the characteristic function of the set {(x,t) £ K" +1 | u > 0} (see Section 
3.4 for the precise definition), & is a given open set in K n+1 , and u is a locally 
bounded weak (i.e., in the distributional sense) solution. 

Observe that {u > 0} is a priori an unknown open subset of 9. We denote by F 
the intersection of 9 with the boundary of the set {u = 0}. We will call F the free 
boundary. For our dimple forming, we have F = {F(f ), {f}},>o of Fig. [5] 

If, additionally, the condition 

^" ^ a.e. in 0, (16) 
at 

is satisfied, then our FBP ( [T5| becomes the Stefan problem mentioned above. In- 



equality ( 16 1 means that our dimple is formed without "returning back". In general, 



however, we can not guarantee ( 16 1 from the assumptions given in Sections 3.2 and 
[331only. 

For the function f(x, t), we assume that 

(1) /is non-degenerate in 9, i.e., there exists 5q > such that f(x, t) > 5o for any 

0,0 s 9; 

(2) / is Holder continuous in 9 with some a G (0, 1), i.e., / is a bounded contin- 
uous function in 2>, and for all points (x,/) and (y,s) € 9 such that (x,t) ^ (y,s) 
we have the inequality 

\f(x,t)-f(y,s)\ 



-3f + |r-s| 2 ) a/2 



4.2 Qualitative Properties of Solutions 



Let m be a solution of FBP ( 15 i, let z° = (x°,t°) € ^, let p > sufficiently small, 
and let Q p (z°) := — x°| <p} x (f° — p 2 ,t°). The following estimates provide 
us with information about the behaviour of our displacement u near the interface 
between the sets {u > 0} and {u = 0}, i.e., near the free boundary. 

• There exists a constant c = c(n) > such that 

sup u ^ c(n)p 2 . (17) 

Q P (z°) 

This nondegeneracy estimate holds true for all points z° belonging to the closure 
of the set {u > 0} and for all p sufficiently small. 

Moreover, our solution u has quadratic growth near the free boundary: 

• There exists a constant C > such that 

sup u ^ Cp 2 . (18) 

Gp(z°) 
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This nondegeneracy estimate holds true for all free boundary points z° € F and 
for all p sufficiently small. 



There exists a universal constant M > such that 



sup 

Gp(z )n{«>0} 



d 2 u 


+ 


du 




dxidxj 






) 



CM. 



(19) 



This inequality holds true for all free boundary points z° G F and for all p suffi- 
ciently small. 

For f(x,t) = const all these three statements were proved in ll47l . The case of 
general / = f(x, t) was considered in [48] for n > 1 and in [49 1 for n = 1. 

The preceding estimates \YT) , ( |18) , and (|T9j indicate that characteristic base di- 
ameter, growth rate, depth and time of the dimple forming are mathematically well 
defined and, therefore, these values should, in principle, be measurable in experi- 
ment, see below Sect. 15.31 



4.3 Classification of Blow-up Limits in W +1 

The idea is to use blow-up sequences, which are a kind of zooms, and to look at the 



"infinite zoom". Suppose that u is a solution of the problem ( 15 1, z* = (x* ,t*) is a 



free boundary point and f{x* , t * ) ^ 0. For X > consider the functions 



u | x" +x—, X J*+t- l ~ 



V/(*v*) J 1 

u x (x,t):=^ jj- for (x,t)e9 x :=j-9. 

There exists a sub-sequence {A^} converging to zero such that the blow-up sequence 
{ux k } converges to one of the following blow-up limits: 

• uo — UQ te {x,t) := j(x r • e) 2 , for a unit vector e, where x T ■ e denotes the scalar 
product in W, 

• uo = uo,m(x,t) = mt +x T ■ ^# • x, where mis a constant and ^ is a (n x n)- 
matrix satisfying Tr^# = m + 1 . 

Observe that the blow-up limits can (in general) depend on the choice of the sub- 
sequence {Aa}. But it should be emphasized that in view of the non-negativity of u, 
the limit function uq is the unique non-negative distributional solution of 

Au-^-=Xu>o a.e. in E" x (-°°,t*). 
dt 

This means that in the second case m and ^# are defined uniquely as well. 

For f(x,t) — const these results were obtained in [47 1. For the general case / = 
f(x, t) we refer the reader to [48 . 
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4.4 Classification of the Free Boundary Points 

Going to the regularity properties of the free boundary we observe that 

• The free boundary is a closed set of zero (n + 1)-Lebesgue measure. 

• The free boundary F = {"regular points"} U {"singular points"}. 

• The set of the singular free boundary points is closed. 

The singular points are defined as the free boundary points for which there ex- 
ists a blow-up limit of the second type, i.e., u°(x,t) = mt + x T ■ ^ ■ x. The set 
r \ {"singular points"} is the set of regular points. For singular points and for 
k = 0, ...,n the sets S(k) are considered additionally, where S(k) is defined as the 
set of singular points such that dim Kern ^ = k and the smallest of the k non-zero 
eigenvalues is bounded from below by a fixed positive constant. 

The complete classification of all free boundary points can be given via a rela- 
tively new approach introduced by G.S. Weiss in |50|. For a solution u of FBP ( 15 i 



and for a free boundary point z* = (x*,t*) consider the following energy functional 

W(x,z*, u,f):= i J j \\Vu\ z +2fu + -^— i J G(x-x*,t* -t)dxdt, 

r*-4-r 2 \x— X*|<T 

where G(x,t) is the heat kernel. 

The functional W has the following remarkable properties: 

• lim W(t,z* ,u,f) exists and is finite; 

T->0+ 

• There are only two possible values for lim W(T,z*,u,f), namely 

T->0+ 

(A n , if z* is a regular point 
lim W (t,z ,u,f) = < 

t^o+ I 2A„, if z is a singular point 

for some constant A„ > depending only on the dimension n. 
Finally, we observe that 

• Around regular points the free boundary is a smooth graph. 

• Singular points belonging to S(n) are isolated. 

• S(n) is contained locally in a C^-graph in space. 

1 12 

• S(k) for ^ k ^ n — 1 is contained locally in a fc-manifold of class C x t . 

For all results concerning the regular points we refer to flTl . The results about 
singular points were proved in [48 1 (see also |3T)). 
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Our conclusions consist of preliminary findings which will require further experi- 
ments and measurements to be confirmed - or falsified. 



5.1 Summary of (Partly Speculative) Working Hypotheses 

For better reading, we begin this closing section by summarizing our assumptions 
and choices of emphasis. 

1 . Among all relevant aspects of the well- and malfunctioning of pancreatic j3 -cells, 
we focus on a single membrane process, the lipid bilayer fusion event. 

2. We suppose that future imaging will prove the making of a plasma membrane 
dimple before the making of the fusion pore. 

3. We assume that the dimple making is essential for the performance of regulated 
exocytosis also in j3 -cells. 

4. We suspect that maintaining the fusion pore and continuing release of the content 
of the insulin granules will be interrupted when the dimple is not preserved. 

5. We argue for a long-distance regulation of the dimple making (and, hence, the 
secretion process), i.e., we claim that the docking of readily releasable insulin 
granules at the plasma membrane and the consecutive making of the fusion pore 
and the release of the hormone molecules is induced not as a purely local phe- 
nomenon and hence can not be explained solely by glucose stimuli and corre- 
sponding Ca 2+ influx through ion channels in the neighborhood of the release 
site. Instead of that, we claim that stimuli and ion influx also far from the exo- 
cytosis site and mediated by intra-cellular signaling and energy transport will be 
decisive for initiating and maintaining regulated exocytosis even at one isolated 
site. 

6. We point to a low-frequent electromagnetic field wave as a possible regulator of 
exocytosis. 

7. We suppose that the electromagnetic field will be closed over the plasma mem- 
brane due to the iron content of enzymes embedded between the phospholipids 
of the plasma membrane. 

8. We assume that the electrical activity of the mitochondria is equally important as 
the electrical activity at the plasma membrane. We suppose that the mitochondria 
sequester and release Ca 2+ ions in a self -regulated way which builds an electro- 
magnetic field and generates a directed (travelling) field wave. 

9. We suppose that the synchronization of the electrical activity of neighboring mi- 
tochondria is due to energy efficiency. 

10. While many aspects of j3-cell function will require an analysis of their collec- 
tive functioning in the Langerhans islets, we conjecture that essential aspects of 
regulated exocytosis can be observed on the level of a single cell. 
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We have provided a mathematical model for the initiation of regulated exocytosis 
and the making of the fusion pore. The model relates the geometry and the dy- 
namics of one single membrane process, namely the forming of an inward oriented 
dimple in the plasma membrane before the fusion event, with electromagnetic fea- 
tures of intracellular calcium oscillations. The model suggests a new explanation for 
the observed flickering of regulated exocytosis, the vanishing of the first phase of 
secretion in stressed or tired j3 -cells and the final halt of all secretion in overworked 
dysfunctional cells: the electromagnetic free boundary model points to the lack of 
stability and coordination of the intracellular Ca 2+ oscillations prior to the bilayer 
membrane vesicle fusion. We recall that the field character of these oscillations is 
magnetic (therefore transferring energy to the fusion site at the plasma membrane 
without any loss). It must be distinguished from the widely studied Ca 2+ influx 
changing the electrostatic potential across the plasma membrane and accompanying 
regulated exocytosis. 

The model is based solely on physical First Principles. All parameters have a 
biophysical meaning and can, in principle, be measured. 



5.3 Suggested Experiments and Measurements 

At the present stage of our knowledge, the mathematical and biophysical correct- 
ness of our model does not prove its relevance for explaining the phenomena it 
claims to explain. To decide whether the here described phenomena and effects are 
dispensable or decisive for regulated exocytosis, the scales of the ion oscillations, 
the electromagnetic fields, the acting forces, the entering material constants and the 
characteristic times and lengths must be determined. 

Hence, the framework of our electromagnetic free boundary problem for the dim- 
ple making suggest the following array of experiments and measurements: 

1 . We shall observe the Ca 2+ oscillations prior to the bilayer membrane vesicle fu- 
sion also in pancreatic j3 -cells and determine their spatial and temporal character. 
In particular, the observations must 

• check the intracellular origin of the oscillations; we may, e.g., deliberately 
silence (empty) some types of organelles by adding suitable agents, see, e.g., 
Fridlyand et al. Il52l : 

• identify the participating organelles (Ca 2+ storages); and 

• decide about the orientation (the direction) of oscillations, the frequency, de- 
pending on stimulus, and the distinction between almost simultaneous oscil- 
lations pointing in different directions as their selected sites for the making of 
the fusion pore. 
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Experiment and Discovery 

In science, we meet various ways of relating experiments and discovery and very different 
perceptions of that process. Physics has (naturally had to) developed what other scientists consider 
a rather crude iteration cycle between theory, model and experiment: first we mess around making 
some observations; then we specify and embed the observation into a mathematical model; after 
which we create a theory by relating our observations and our model to other observations and 
other models; all that gives rise to new questions and new experiments, partly model-based, partly 
again messing around - and the cycle begins again. 

A related view, shared in systems biology, upon experiments is to distinguish between 
exploratory experiments, anteceding models and theory, and critical experiments confirming or 
falsifying a model or a theory. For use in medicine and biology, the immunologist and 1960 Nobel 
Laureate Peter Medawar elaborates on that distinction and discerns four kinds of experiments. 

A first kind of experimentation is one in the original Baconian sense, i.e., any organized 
experience or happening, devised with ingenuity or skill. It is the consequence of "trying things 
out" or even of merely messing around. Baconian experiments try to answer the question "I wonder 
what would happen if ..." Such experimentation proved productive from the time of alchemy until 
the discovery of high temperature superconducting a few decades ago. 

A second kind of experimentation is what Medawar baptized Aristotelian: these experiments are 
also contrived, but not out of curiosity but for supporting an existing model or theory or just to 
demonstrate its basic essence by an observatory example, like Galileo's observation of the Jupiter 
moons in confirmation of the Copernican planetary system and NASA's recent Lunar Laser 
Ranging Test of the invariance of c. 

The ideal of modern science is a third kind of experimentation, the Galilean or critical 
experiment. The main goal is to provide striking evidence against a previous erroneous perception. 
The aim is clarification by falsification, like Galileo's balls rolling down an inclined plane in an 
almost frictionless environment to refute mass dependence of the velocity of falling bodies (and 
thereby establishing constant acceleration) or the famous Michelson-Morley experiment to reject 
the theory of a luminiferous aether. 

Finally, there are the Kantian experiments - the thought experiments: "Let's see what would 
follow if we took a somewhat different view..." Kantian experimentation requires no apparatus 
except sometimes a computer. From discussions between Bohr and Einstein, series of Kantian 
experiments are documented. Einstein himself attributed his discovery of Special Relativity to 
shared fantasy with his first wife of sitting together on the tip of a light beam. 

More to read: P.B. Medawar, Advice to a Young Scientist, Basic Books Inc., New York, 1979. 

Added by the editors 



2. We shall modulate the oscillations by submitting the cells to an external field 
generator with variable frequencies, to prove the magnetic character of the field 
wave associated to the oscillations. 

3. We shall determine the surface tension in the plasma membrane of living cells. In 
particular, we shall measure and/or calculate the bending rigidity and stretching 
elasticity under "repacking" of the ball shaped heads of the inner lipids under 
area changes. For living cells, we expect that these magnitudes are substantially 
larger than for model membranes, e.g., due to osmotic pressure in living cells, 
see Henriksen and Ipsen (53]. 

4. We shall measure the cytosol viscosity close to the plasma membrane, i.e., update 
the classic study [54] by Bicknese et al., and locate actin filaments blocking for 
unwanted docking of the insulin granules at the plasma membrane. 
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5. We shall check Rosenheck's estimate (I.e.) for charged molecules in the plasma 
membrane. 

6. We shall estimate the content of (para-, not necessarily ferro-)magnetizable Fe 
atoms and crystals in the membranes to determine their magnetic momentum. 

7. We shall estimate the distribution of the inhomogeneities of the magnetic field 
near the plasma membrane. 

8. We need precise electron or atomic force microscope slices of the dimple making 
and of the degree of its singularity in j3 -cells. 

9. We shall measure by patch clamp technique the expected decrease of the capaci- 
tive reactance X c under dimple forming. 

10. We shall correlate the Ca 2+ oscillations with the fusion events; in particular, we 
shall confirm the spatial and temporal coincidence of flickering of exocytosis 
with break-downs of the field wave. 

11. We shall demonstrate the absence or weakness of the Ca 2+ oscillations after 
stimulation in stressed or tired j3 -cells. 

12. After obtaining reliable values of all data involved in our mathematical model, we 
shall create a computer simulation of the free boundary problem to get a solution 
graphically (in form of a surface). After that, we can compare the dimple images 
and the modeled surface. 
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